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ABSTRACT 

A new method is proposed for modelling spherically symmetric inhomogeneities in the 
Universe. The inhomogeneities have finite size and are compensated, so they do not 
exert any measurable gravitational force beyond their boundary. The region exterior to 
the perturbation is represented by a Friedmann-Robertson- Walker (FRW) Universe, 
which we use to study the anisotropy in the cosmic microwave background (CMB) 
induced by the cluster. All calculations are performed in a single, global coordinate 
system, with nonlinear gravitational effects fully incorporated. An advantage of the 
gauge choices employed here is that the resultant equations are essentially Newtonian 
in form. Examination of the problem of specifying initial data shows that the new 
model presented here has many advantages over 'Swiss cheese' and other models. 
Numerical implementation of the equations derived here is described in a subsequent 
paper. 

Key words: Gravitation - cosmology: theory - cosmology : gravitational lensing - 
cosmic microwave background - galaxies: clustering 



1 INTRODUCTION 

Many attempts have been made to model secondary 
anisotropies in the cosmic microwave background (CMB) 
due to galaxy clusters. The most widely employed approach 
is to us 
tional 
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regime; and all observations can be modelled in the exterior 
region, providing clear, unambiguous predictions of the ef- 
fect of the cluster. The main disadvantage of such models is 
that the matter distribution is unrea listic, and appears to 
overe stimate the effect of the cluster (Quilis, Ibanez & Saez 
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199$ . 

More realistic models can be constructed by working 
with an arbitrary density profile while restricting to spher- 
ical symme t ry and ignoring pressure ( |Panek 1992 ; |Arnavj 
et al. 199J;_feaez, Arnau & Fullana 1993; Arnau, Fullana 



Saez 1994 



Fullana, Saez & Arnau 1994; Quilis & Saez 



199S ) . Spherically symmetric models with vanishing pres- 
sure have implicit analytic solutions in General Relativ- 
ity, one simple f orm of which is provided by the Tolman- 



Bondi solution (Iblman 1934; Bondi 1947). Models based 



full, nonlinear theory employed 'Swiss cheese' models ( 


Reef 


& Sciama 1968; Dyer 1976; Kaiser 19821 


Nottale 1982a 




Not- 


tale 1982b; Nottale 1983; 


Nottale 1984) 


. In these models a 



collapsing FRW Universe is surrounded by a compensating 
vacuum region, which then matches onto an exterior ex- 
panding FRW Universe. The compensating region ensures 
that the perturbation has no net gravitational effect on the 
exterior Universe. These models have two key advantages: 
analytic calculations can be performed in the fully nonlinear 
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on this solution provide a more reliable method for comput- 
ing anisotropies due to clusters or voids, and it is a hybrid 
of this scheme that we discuss here. Previous work with 
Tolman-Bondi models has considered density and velocity 
perturbations, but compensating effects are not usually in- 
cluded. (One excepti on to this is the 'type F class of models 
considered by Panek (1992), but these are only compensated 
at infinity.) The fact that observers comoving with the fluid 
are not in a region modelled by a homogeneous FRW cos- 
mology makes it harder to discuss perturbations, since these 
observers see a dipole anisotropy in the CMB, which is es- 
sentially an artefact of the model. 

One difficulty in constructing compensated models is 
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that the initial density and velocity profiles must be chosen 
in such a way that streamline crossing is avoided. Without 
this, shock fronts form and one would have to include pres- 
sure to produce a realistic model. For example, if one just 
perturbs the density profile without perturbing the velocity 
field, then streamline crossing is inevitable in any compen- 
sated model other than those of 'Swiss cheese' type (where 
the problem would only occur in the vacuum region). 

Here we discuss a family of models which avoid the 
problem of streamline crossing in a very simple manner, 
while keeping the density profile compensated and realistic. 
The initial perturbation is of finite extent and is controlled 
by two physical parameters, one describing the magnitude 
of the perturbation and the other its linear extent. The per- 
turbation drops to zero at a finite radius, beyond which the 
Universe is described by an FRW model. A third param- 
eter, m, controls the degree of the polynomial describing 
the perturbation and can be adjusted to produce the de- 
sired structure in the nonlinear regime. The two parameters 
controlling the perturbation and the further two parame- 
ters describing the external universe uniquely determine the 
form of the velocity perturbation for each polynomial degree. 
This is achieved by matching the first m derivatives at the 
centre and edge of the perturbation, leaving a polynomial of 
degree 2m + 1. From the velocity profile the density pertur- 
bation is uniquely fixed by the constraints that there are no 
decaying modes present and that the density distribution is 
compensated. In these models density compensation holds 
at all later times and the velocity field evolves in a way that 
avoids streamline crossing. The fact that the initial density 
and velocity profiles are described by simple polynomials 
also simplifies the task of carrying out accurate numerical 
work. 

In this paper we describe the theoretical aspects of our 
model and provide a detailed comparison with other work. 
In Section y we introduce a novel way of treating the field 
equations for spherically symmetric dust, which arose from 
a new, gauge-theoretic treatment of gravity presented else- 
where (Lasenby, Doran & Gull 1998). This method works 
with a single global time coordinate which reduces the equa- 
tions to an essentially Newtonian form. The insights pro- 
vided by this gauge lead naturally to the idea of specifying 
initial conditions on a surface of constant 'Newtonian' time, 
and we show how to parameterise the matter streamlines to 
take advantage of this possibility. 

In Section ^ the equations for photon propagation are 
presented, and we derive a simple formula for the tempera- 
ture decrement induced in the CMB by the cluster. We also 
derive a new formula for density evolution along a stream- 
line. This considerably simplifies the task of computing the 
density profile at later times. The advantages of our gauge 
choices are further illustrated by linearising the field equa- 
tions to produce a relation between the velocity and density 
profiles for realistic initial conditions. In Section ^ we discuss 
Swiss Cheese models and give a new formula for the CMB 
temperature decrement in such models. We also consider 
the strengths and limitations of these models as a prelude 
to presenting an improved scheme. 

Our improved model is presented in Section ^. De- 
spite the simplicity of the model, numerical simulations re- 
veal that the initial perturbation evolves to produce a den- 
sity profile which accurately matches the King profile ob- 



served in many galaxy clusters. The main limitations of 
our model concern the imposition of spherical symmetry 
and the lack of pressure support in the evolving cluster. 
Both assumptions lead to large cluster infall velocities rel- 
ative to the CMB, and these are certainly much larger 
than the velocities obs erved in nearby clusters. How ever, 
recent numerical work (Navarro, Frenk & White 1996) has 



suggested that large infall velocities might occur in evolv- 
ing clusters at higher redshifts, for which our model would 
then be well suited. For example, the microwave decre- 
me nts as sociated with quas ar pai rs reported by Jones et 
al. ( 1997 ) and Richards et al. ( 1997 ) may be due to combined 
Sunyaev-Zel'dovich and gravitational effects from high red- 
shift cl usters. This is discusse d further in an accompanying 
paper ( Dabrowski et al. 1997 ), where the theoretical model 
for the gravitational effect described here is implemented nu- 
merically. Some results of this analysis are briefly discussed 
in the concluding section. Natural units G = c = Th = 1 are 
employed throughout. 



2 THE FIELD EQUATIONS AND THE 
NEWTONIAN GAUGE 

The equation s use d here were first derived by Lasenby, Do- 
ran & Gull (1995) using a new, gauge-theoretic approach 
to gravity. This approach employed 'geometric algebra', the 
mathematical language which seems to capture best the na- 
ture of the gauge fields and the structure of the equations. 
Despite the different concept ual fo undations of the theory 
developed by Lasenby et al. (1998), the predictions of the 
theory agree with general relativity for a wide range of phe- 
nomena, and for the systems studied here the two theories 
agree exactly. So as not to alienate readers unfamiliar with 
geometric algebra, we have not employed it in this paper, 
and have instead adopted an approach closer to that of 
standard general relativity. The price of this approach is 
that some of the results em ploye d below are quoted without 
proof from Lasenby et al. (1998). 



The system we study in this paper consists of a spher- 
ically symmetric distribution of dust. Modifications to in- 
clude either pressure or a cosm ologi cal constant are not 
considered here. Lasenby et al. ( 199S ) showed that a nat- 



ural gauge emerges for the analysis of such systems. This 
gauge is global, employing a single time coordinate t and 
a radial coordinate r. The time coordinate t measures the 
free- fall time for observers comoving with the fluid. In inho- 
mogeneous regions the radial coordinate r is related to the 
strength of the tidal force defined by the Riemann tensor. 
The key dynamical variables are the density p(t, r), a veloc- 
ity field u(t, r) and a gene ralise d 'bo ost' fa ctor denoted by 
T(t,r). (In Lasenby et al. ( [l998| ) and ( |l995| ), T was denoted 
by gi and u was denoted by g2.) The fields V and u define 
the line element by 



ds 



r 2 



r 2 



dr 



-r 2 (d6 2 +sm 2 9d(j> 2 ), 



(1) 



where 6 and <f> are standard polar coordinates, < 6 < n, 
< 8 < 2n. The density p(t, r) defines a total gravitational 
mass M(t, r) by 
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M{t,r) 



4ns p(t,s)ds, 



and the variables F, u and M are related by 



r- 



2M/r + u 



(2) 



(3) 



If we rewrite this equation in the form 

iu 2 -Af/r=I(r 2 -l) (4) 

we see that it has a simple, Newtonian interpretation: it is 
a Bernoulli equation for zero pressure and non-relativistic 
energy (r 2 - l)/2. 

The variable u defines the integral curves of the con- 
served fluid current by 

!=«■ * 5 > 

These integral curves are also matter geodesies, since there 
is no pressure present. The key dynamical equations are up- 
date equations along these fluid streamlines, for which we 
require the comoving derivative 

D - 9 , 9 (R\ 

Di = di +U d-r- (6) 
The total gravitational mass enclosed, M, is conserved along 
the fluid streamlines, 



DM 
~Dt 



0. 



(7) 



Implicit here is the assumption that the initial conditions 
are chosen to avoid the possibility of streamline crossing. 
The Euler equation in this setup is 



Du 

~Dt 7" 
It follows that 



M 

72" 



(8) 



= 0. 



(9) 



DT 
~Dt 

The above set of equations look Newtonian, yet they are 
fully consistent with general relativity. This is the advantage 
of the gauge choices made to manipulate the equations into 
the above form. The form of the equations led us to name 
this gauge the 'Newtonian' gauge. Initial data for the above 
set of equations can be given in a variety of forms, but the 
most natural is to specify the density p(ti, r) and the velocity 
u(ti, r) at some initial time t,. From these one can calculate 
M(ti,r) and T(ti,r), which are then conserved along the 
streamlines. 

To complete the set of basic equations it is useful to 
introduce the velocity gradient H: 

H(t,r) = ^-u(t,r). (10) 
or 

From the above equations it is straightforward to show that 
^ = -(2u/r + H)p. (11) 

The function H is shown below to play a fundamental role 
in determining the temperature fluctuations caused by an 
infalling cluster. 

The above set of equation s hav e analytic solu tions which 
were first found by Tolman ( |l934| ) and Bondi fll947| ). The 
difference between the above approach and the Tolman solu- 
tion is that Tolman used M as the radial coordinate, instead 



of r. This has the disadvantage of hiding the Newtonian na- 
ture of the equations, making it hard to visualise the physics. 
The form of the analytic solution to the above equations de- 
pends on the sign of T 2 — 1 on the initial timeslice. There 
are three cases to consider: 



i. r 2 < 1 

For this case the matter streamlines (geodesies) are defined 
by 

M 



1 



t-U 



r 2 

M 



(1 — cos^) 

(■q — sin?? — rji + sinr/i) 



(12) 
(13) 



(i_r 2 )3/2 

where r\ parameterises the curve and r\i is determined from 
the initial value of r at time ti. The quantities M and V arc 
constant along the curve. In order to choose the initial value 
of r\ correctly, one needs the further result that 

M 



r . (1 _ r2) i/2 sil »7 ( 14 ) 

so that T)i is uniquely fixed between < rji < 2n from a 
knowledge of whether the initial velocity is inwards or out- 
wards. The case rji — n cor respo nds to starting from rest, as 
treated by Lasenby et al. (1998). 



2. r 2 = 1 

This case includes a flat FRW cosmology, and the equations 
can be integrated directly to give 



t-u 



2 ( r 3/2 _ r .3/2) 

3 (2M)V2 • 



(15) 



The velocity is chosen to be outwards to avoid a singularity 
forming instantaneously. 

3. T 2 > 1 

This case includes open cosmologies. The matter curves are 
parameterised as follows: 

M 



-(cosh?7 — 1) 



t-U 



M 



(r 2 - 1) 3 / 2 



(sinhr; - 77 - sinh?^ + r]i) 



and the velocity is given by 
M 



r(r 2 - l) 1 / 2 



sinh r\. 



(16) 
(17) 



(18) 



For this case it is also necessary to start with an initial 
outward velocity to avoid streamline crossing. 

The above parameterised expressions for the mat- 
ter stre amlines are more general than those usually pre- 
sented (Panek 1992). This is necessary so that the initial 



data can be specified on a constant t timeslice, which pro- 
vides better control over the form of the initial perturba- 
tions. The initial time U is arbitrary and is usually either 
set to zero or chosen so that the present epoch corresponds 
to t = 0. Initial data is defined by p(ti,r) and u(ti,r) and 
from these the initial values of M and V are found. At the 
initial point r(ti) the value of rji is found, and the resulting 
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curves are plotted parametrically. Examples of t hese c urves 
for collapsing dust are given by Lasenby et al. (1998). The 



lack of pressure support means that infalling matter col- 
lapses to a singularity in a finite time. This does not pose 
any problems, provided photons passing through the cluster 
do not encounter a singularity or horizon. 



3 PHYSICS IN THE NEWTONIAN GAUGE 

There seems little doubt that the Newtonian gauge provides 
the cleanest description of gravitational phy s ics for spheri- 



cally s ymmetric systems (Lasenby et al. 1998; Lasenby et al 



1997a) and it is surprising that it is not more widely used. 



Here we review here some key physical results from the New- 
tonian gauge perspective. These results are employed in later 
sections to construct numerical models. A further illustra- 
tion of the use of the Newtonian gauge which is slightly 
outside the main theme of this paper is contained in Ap- 
pendix [a|. 

3 . 1 Cosmology 

While the Newtonian gauge arises naturally in the study of 
spherically symmetric systems, it is not immediately appar- 
ent that it is valuable in the study of homogeneous cosmo- 
logical models. Since the models we discuss below contain a 
homogeneous region outside the perturbed region, it is useful 
to see how standard cosmological notions fit into th is sch eme 
(see also Gautreau ( 1984 ) and Ellis & Rothman ( 1993 ) for 
similar, more detailed discussions). 

In a homogeneous model the density p is a function of 
time t only, so we have 



M(r,t) = -7rr 3 pW- 



(19) 



Equations (j7|) and ( |ll| ) now yield H = u/r and p — —3Hp, 
where the overdot denotes differentiation with respect to 
time. It follows that H is a function of time only for homo- 
geneous models. 

The Euler equation (|^) yields the second cosmological 
equation, 



H + H 2 = -fp. 



(20) 



It is clear now that we can identify H as the Hubble function 
and that the Newtonian time parameter agrees with cosmic 
time in homogeneous regions of the model. 

To complete the set of cosmological equations we note 
from equation (tfl) that T 2 must be of the form 



V 2 = l + T 2 X {t). 



(21) 



Conservation of V down the streamlines (equation (|9|)) shows 
that 



X = -2H(t) X , 
hence we have 



r 2 = 1 - fcr 2 exp{-2y H(t')dt'}. 



(22) 



(23) 



The Friedmann equations in their standard form (with zero 
pressure and cosmological constant) are recovered by intro- 
ducing the scale factor S(t) via the definition 



H{t) 



s(ty 



so that 



T 2 = 1 - kr 2 /S 2 , 



(24) 



(25) 



with k = 0, ±1 determining the type of cosmology in the 
usual manner. 

In the Newtonian gauge version of cosmology, particles 
in the Hubble flow are pictured as moving radially outwards 
at a velocity r = H(t)r. The expansion centre is not a phys- 
ical feature, since the model is homogeneous and all gauge- 
invariant (measurable) quantities are functions of time only. 
One can often take advantage of this fact by considering ob- 
servers located at the origin. This provides a novel way of 
visualising some of the less in tuitive aspects of FRW cos- 
mologies ( Lasenby et al. 1998 ). Along the matter stream- 
lines r(t)/S(t) is constant. The standard 'co-moving' gauge 
for cosmology is obtained by introducing a new radial co- 
ordinate R — r / S, so that the matter geodesies are lines of 
constant R. 



3.2 Photon Paths and Redshifts 

In studying photon trajectories it is sufficient to just con- 
sider motion in the 9 = 7r/2 plane. We can parameterise the 
photon paths in terms of the Newtonian time parameter t 
(cosmic time for homogeneous models), so that the path is 
defined by r(t) and 4>{i). The condition that the trajectory 
is null means that we can write 



dr 
~dt 



r cosx + u 



(26) 
(27) 



dip smx 
~dt ~ T 

where x is the angle between the photon trajectory and the 
cluster centre, as measured by observers comoving with the 
fluid. The geodesic equation is equivalent to the following 
first-order equation for x- 

^ = sinx(~r/r + (H — u/r) cos X ) . (28) 

The above set of three first-order equations is simpler to 
implement numerically than the second-order equations ob- 
tained from the direct approach. The equations have the 
further advantage of dealing directly with the measurable 
quantities t and x- It is typical of the gauge theory approach 
that one obtains first-order equations in the gauge invari- 
ant observables such as these. Initial data for the trajectory 
consists of position data n and <pi and a direction Xj- For 
most calculations, however, the data is given in the form of 
the observer's position and an angle on the sky x, and the 
equations are then run backwards in time to take the pho- 
ton back through the cluster. This is the simplest way to 
perform lensing simulations. 

The remaining content of the geodesic equation con- 
cerns the frequency to in the rest frame o f the fluid. This 
satisfies the equation (Lasenby et al. 1998) 



du> 
~dt 



~uj(H cos 2 x + — sin 2 x) 



(29) 



It is not hard to show that the above equations imply that 
the angular momentum L = u>r 2 4> = — ruismx is conserved. 
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For the case of a homogeneous cosmology, the redshift 
equation rt29|) reduces to 



dui 



-uH, 



(30) 



from which we recover the standard result that ojS is con- 
stant for a photon travelling in a homogeneous, dust-filled 
Universe. This result shows how the standard predictions of 
cosmology are recovered in this unfamiliar gauge. 

To calculate the effect of a cluster on the CMB, we write 



u(t,r) = rH e (t) + A{t, 



(31) 



where H e (t) is the Hubble function in the exterior Universe 
at the same time t, and A is the difference between the equiv- 
alent velocity of the unperturbed Universe and the cluster 
velocity. For a photon passing through the cluster we find 
that 



d 



ln(wS) = 



dA 2 
■ cos x 



A . 2 

-sin X- 



dt"' x ~~' dr 
We therefore define the function 



T 2 , ^A 

J il 



cos 2 x+ — sin 2 x), 



(32) 



(33) 



where the integral is evaluated along the photon path be- 
tween the time the photon enters the cluster (ti) and the 
time it leaves {t2)- The function e is small, since the contri- 
bution to the integral from near the cluster centre tends to 
cancel the contributions from further out. The main effect 
producing a non-zero e is essentially the evolution of A with 
time. 

If the photon enters the cluster from the CMB with 
frequency 101 and reemerges with frequency ui'2 we have 



S(ti) - 



(34) 



An equivalent unperturbed photon will have frequency 
S(ti) 



U>2 = OJl 



Site) ■ 



so the physically measurable temperature decrement is 

AT 
~~f~ 



U) 2 — UJ2 



1 



(35) 



(36) 



This simple result enables easy computation of the effect of 
the cluster for various angles on the sky. 



3.3 Evolution of the Physical Variables 

The two key physical variables required at later times are 
the density distribution, which should evolve to produce a 
realistic distribution for a cluster, and the velocity gradi- 
ent, which controls the propagation of photons through the 
cluster. One way to calculate these at later times is to differ- 
entiate numerically the mass M(t, r) and velocity, which are 
easily found on any given streamline. A better approach, 
however, is to construct explicit analytical formulae for p 
and H on a given streamline, which we now discuss. 

We start with the density profile and employ the rela- 
tion 



p(t ' r) = (^) t ^ P(ri) ' 
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(37) 



which gives p(t, r) as a function of the initial density profile 
piji) = p(ti,ri) and the derivative of the streamline starting 
position r;(f,r) with respect to r. To calculate this latter 
term we use the reciprocity relation 



(drj\ (dr\ I dt \ 
\dr)t \dt) ri KdrJ 



-1 



to deduce that 



\arJt \oriJ 



(38) 



(39) 



The final term can be calculated directly from the param- 
eterised form of the streamlines. It is convenient here to 
introduce the function 



f( ri ) = r(u, n ) 2 -1. 
With this we find that, for / 
1 



(40) 



(-Y 1 



2f(n)M(r 
-87rn 2 p{ri)f(n)(r - ut) 

uM(n) d ( nf(n 



dri 



f(n)u{n) dn \ M(n) 



(41) 



where M(ri) = M(ti,ri) is the initial mass function. With 
this result one need only compute df/dri from the initial 
data to find the value of p on a given streamline. The case 
of / = must be treated separately and yields simply 



(— V 

KdrJt 



1/2 



ut 



2M(r 



-4nn p{ri). 



(42) 



The same technique is used to find the velocity gradi- 
ent on a streamline at later times. On differentiating equa- 
tion (H) we find that 



du 



M 



+ 



\drj t 



l dM(n) | ldf(n 
r drt 2 dri 



(43) 



which again gives if on a streamline in terms of the stream- 
line label ri and the initial data. These formulae are very 
useful for carrying out accurate numerical simulations. 



3.4 The Linearised Equations 

We now derive a physical constraint on the initial conditions 
for our model by considering evolution in the linear regime. 
Linearising the field equations around a homogeneous back- 
ground cosmology is a straightforward exercise in the New- 
tonian gauge. We start by introducing the variables 



Su = u — rH{i) 
Sp = p- p(t) 



SM 



Aits Sp(t, s) ds, 



(44) 
(45) 

(46) 



where the barred quantities denote the background homoge- 
neous values. The linearised equations for the velocity field 
now yield 



DSu 
Dt 
and 



SM 



H8u 



(47) 
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D (D5u\ 
Dt V Dt ) 



+ 3H 



DSu 
Dt 



H i Su = 0, 



(48) 



where we have assumed that the background is described by 
a spatially flat, A = cosmology. Equation ( fjjj ) shows that 
the evolution of the density perturbation down the stream- 
line is determined entirely by the background cosmology, 
and not by the local properties of the perturbation. It fol- 
lows that there is nothing special about the assumption of 
spherical symmetry and equation (j48|) is a general result for 
the growth of a velocity perturbation in a dust-filled cosmol- 
ogy- 

We can parameterise Su in terms of a streamline start- 
ing point and cosmic time t. In this way the streamline 
derivatives are replaced by derivatives with respect to t, and 
we arrive at the equation 



Sil + 3H8u~ H Su = 0. 



(49) 



Since H — 2/(3t) we see that Su has modes going as t l/3 and 
£-4/3 Only the former of these is growing with respect to 
the background value. If we assume that we are in an epoch 
when the decaying mode can be ignored, we find that 

It follows that the density and velocity perturbations in the 
linear regime, when no decaying modes are present, are re- 
lated by 



-HSi 



SM 



(51) 



This provides a constraint on the initial conditions for our 
models. 



4 SWISS CHEESE MODELS 

The most popular method to date for theoretical calcula- 
tions of CMB anisotropics in the nonlinear regime has been 



to employ 'Swiss Cheese' models (Flees & Sciama 1968 



Not- 



tale 1982a). In these models an overdense collapsing region 
is surrounded by a compensating vacuum region which then 
matches onto the external Universe. Nottale (1982a) de- 



veloped a complicated set of equations for describing pho- 
ton redshifts in such a model. The interior region is mod- 
elled by a collapsing fi > 1 cosmology, and three coor- 
dinate systems are employed with two sets of matching 
equations at the boundaries. Initial data consists of the set 
{Ho, po, H c , p c ,ri, z c }, where Ho and po are the values of the 
Hubble constant and density at the present epoch, H c and 
p c are the equivalent variables for the central cluster region 
as the photon passes through the centre, z c is the redshift of 
a photon emitted from the centre of the cluster at the epoch 
described by H c and p c , and is the cluster radius when 
the photon exits the central core. One of the five dimen- 
sionful quantities is arbitrary and can be scaled to unity, so 
Nottale's setup defines a five-parameter model. To analyse 
this it is convenient to introduce a set of dimensionless vari- 
ables. Along with z c these consist of the two deceleration 
parameters q c and qo, together with 



y = H c ri and h = — -. 

Ho 

It is also convenient to define 



(52) 



1/3 



1/3 



,2/3 



(53) 



From the initial data it is possible to calculate the red- 
shift of a photon as it reemerges into the exterior Universe 
given its redshift at the centre z c . The method relies on the 
fact that the initial data define the epoch when the photon 
is at the cluster centre. First H c , p c , r» and z c are used to 
calculate the redshift of the photon at the edge of the col- 
lapsing region. This s tep is a straig htforward application of 
the Mattig relation ( Peebles 1993 , Chapter 13). Next one 
computes the photon redshift after it crosses the void region 
and enters the exterior Universe. This calculation is possible 
because of the existence of a timelike Killing vector in the 
vacuum region, which ensures that oj(F + it) is constant as 
the photon crosses the vacuum. From this it is possible to 
calculate the redshift of the photon as it reenters the un- 
perturbed Universe as a function of the input variables. A 
similar, more involved calculation is performed in the oppo- 
site time direction to obtain the total photon redshift across 
the perturbed region. This redshift is then compared with 
that of a photon coming from an unperturbed region of the 
Universe to obtain a central temperature decrement. 

The resultant equations are com plicated and cannot be 
solved analytically. Nottale (1982a) obtained a power se- 



ries solution by expanding in both y (which is necessarily 
small) and z c . A better approach, however, is to obtain an 
expression valid for general z c using only an expansion in y. 
This calculation is tricky, but with the help of the symbolic 
algebra package Maple it is possible to perform such an ex- 
pansion. The result, which we give here for the first time, 
is 



~F = ~Mo(^fY {(l + z c )(l + 2q z c ) 1 

^-'"(ri 1 -)]-"}' 



(54) 



That one can derive an equation such as ( |54[ ) is some- 
what surprising. From the Newtonian gauge perspective, the 
initial data appears to under-determine the system. This is 
illustrated in Figure |l[ which shows that the velocity profile 
is not uniquely specified since various curves can be fitted in 
the vacuum region. Different curves produce different cross- 
ing times for photons to cross the vacuum region and so 
produce different physical models. The trick which makes 
the derivation of ( |54[ ) possible is that z c is specified as part 
of the initial data. This effectively places an extra constraint 
on the form of the velocity function in the vacuum and, while 
not pinning down the form of the curve uniquely, this does 
render the predictions unique. The residual freedom in the 
choice of velocity function amounts to a gauge freedom so 
has no physical significance. 

While the result ( |54| ) provides a useful estimate of the 
effect of a cluster, this type of Swiss Cheese model has a 
number of weaknesses. The density distribution is highly 
unrealistic and must limit the accuracy of any predictions 
from such models. The same can be said of the the high 
infall velocities in such models, which appear to overesti- 
mate the CMB temperature decrement. A further weakness 
is the number of arbitrary parameters in the model. While 
the cluster width, density and redshift are certainly sensible 
parameters for such a model, the velocity in the form of H c 
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5 SIMPLE, FOUR-PARAMETER MODELS 
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Figure 1. Initial Data for Swiss Cheese Models, (a) shows the 
initial density distribution, with a uniform central overdensity 
and a compensating vacuum. The width of the compensating re- 
gion ensures that the total mass enclosed is the same as would 
have been enclosed if the density was uniformly distributed with 
the exterior value, (b) shows the initial velocity distribution. The 
central region is modelled by a collapsing Q > 1 Universe, and 
the exterior region by an arbitrary, expanding model. There is no 
unique curve extrapolating the initial velocity distribution across 
the vacuum. Two different curves are shown, which would pro- 
duce two different models. 



or h is not an obvious physical parameter. No ttale addressed 
this weakness in a later paper (Nottale 1984) by introducing 



the additional restriction that the two matter regions issued 
from the same singularity. This correlates the (Newtonian) 
times in the internal and external regions, and allows one 
to determine h in terms of the remaining dimensionless pa- 
rameters. This assumption is partially supported by Silk's 
result that the effect of eliminating the decaying modes is to 
set the time that each s treamline leaves the initial singular- 
ity to a constant value (Silk 1977). This result only holds in 



the linear theory, however, so it is not obviously applicable 
to Swiss Cheese models. In the following section we propose 
an alternative model which addresses these weaknesses in a 
simple and elegant manner. 



The control over the initial conditions afforded by the New- 
tonian gauge approach provides for many models of the de- 
velopment of a perturbation away from a homogeneous cos- 
mology. Here we discuss a family of simple four-parameter 
models based on polynomial perturbations in the density 
and velocity fields. The perturbation is of finite extent and 
the density distribution is compensated. The region outside 
the perturbation therefore evolves as a homogeneous cos- 
mology, and placing observers in this region allows for un- 
ambiguous calculations of the CMB perturbation caused by 
the cluster. The perturbation is assumed to have grown from 
primordial fluctuations in the very early universe, and to 
still be well described by the linearised Einstein equations. 
We therefore expect that the decaying modes are negligible 
at the epoch when we specify the initial conditions. This 
physical requirement is imposed by ensuring that the initial 
velocity and density profiles satisfy equation (fsij). 

In our model the velocity perturbation is controlled by 
two physical parameters — one specifying the width of the 
perturbation, and the other the velocity gradient at the ori- 
gin. For models of cluster formation the velocity gradient at 
the origin is slightly less than that of the unperturbed Uni- 
verse. Voids can be modelled with a slightly greater velocity 
gradient. A third parameter m controls the degree of the 
polynomial describing the perturbation, and is chosen to to 
ensure a realistic density profile in the nonlinear region. The 
polynomial describing the velocity perturbation has degree 
2m + 1 and is formed as follows. At the centre we have u = 
and the velocity gradient is given. All other derivatives up 
to order m are set to zero. At the edge we match u and 
its first m derivatives to the exterior value, u = rHi. These 
boundary conditions specify a unique curve for each value 
of m. Curves with m—2 and 4 are shown in Figure ^. 

With the velocity profile determined, the density profile 
is found from equation (|5l]) . The initial data for the exterior 
universe consists of a value for the Hubble function Hi and 
the density pi. With the perturbation width given by r\, we 
can introduce the dimensionless variables 



v{x) 



u(r,U) - rHi 



P(r, U) 



(55) 



TiHi 

In terms of these, equation ( |5l] ) takes the form 

x 2 f(%) = ~-l-( x2v ( x ))- (56) 
dx 

Since v(x) is a degree 2m + 1 polynomial, f(x) is a poly- 
nomial of degree 2m. Computing the density distribution in 
this manner ensures that it is compensated at the bound- 
ary. The external values of pi and Hi, together with the 
magnitude and width of the perturbation, form a set of four 
parameters, one of which is arbitrary because of the freedom 
to rescale. The initial density and velocity perturbations are 
given by polynomials, so the function f(rt) — T(r^) 2 — 1 also 
has a simple polynomial expression. It follows that all of the 
terms on the right-hand side of equation ( fll| ) can be cal- 
culated algebraically, enabling accurate computation of the 
density on a given streamline. 

An example of the streamlines resulting from such an 
initial setup is shown in Figure ^| The perturbed region has 
a slight undervelocity, so after a finite time it starts to rec- 
ollapse and form a compact cluster. The figure illustrates 
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Figure 2. Initial data for the 4-parameter model. The solid 
curves are for m = 4 and the dashed curves for m = 2. The 
perturbation has radius n = 1 (arbitrary units), beyond which 
the velocity u(ti,r) is proportional to r and the density is uni- 
form. The perturbation is defined by the velocity gradient at the 
origin (set to 0.6 here) and the radius r;. The density profile is 
compensated, so the total mass enclosed in r < r, is the same as 
would have been enclosed if the initial density was uniform. 



that streamline crossing does not occur in our models. An 
example of the density distribution resulting from growing 
this perturbation is shown in Figure ^| This density distri- 
bution cl osely matches the Kin g profile observed in galaxy 
clusters (Dabrowski et al. 1997), a remarkable result given 



the simplicity of the model. 

The fact that the density distribution evolves to a re- 
alistic profile makes it possible to choose input parameters 
to achieve realistic simulations of the effect of a cluster. A 
photon can then be sent through the cluster to simulate the 



Figure 3. Matter streamlines for an m = 3 model. The pertur- 
bation has initial width 1, with Hi = 1 and pi = 3/(87r). The 
velocity gradient at the centre of the perturbation is 0.95. The 
central region is moving inwards relative to the Hubble flow, so 
recollapses to a singularity in a finite time. The form of the initial 
data ensures that streamline crossing does not occur. 



effect on the CMB. The epoch at which the photon enters 
the cluster region is the fourth physical parameter in the 
model. One can then perform detailed calculations of the 
lensing properties of the cluster and its effect on the CMB 



These are c ontained in an accompanying paper (Dabrowski 



The main limitation with the model presented here is 
the lack of pressure support. This has two effects. First, 
the infall velocities in the cluster are unrealistically large 
compared to a virialised cluster. This will tend to overesti- 
mate the temperature decrement due to the cluster, since it 
is the evolution of the velocity difference A(t, r) with time 
which is central in determining the magnitude of the pertur- 
bation. The second effect is that a central singularity forms 
too soon, when the first streamline reaches the origin. There 
are no difficulties in extending the model beyond this time, 
provided the photon has exited the cluster before the hori- 
zon has formed, but the initial conditions must ensure that 
singularities do not form at early epochs. It turns out that in- 
troducing the perturbation at recombination does not pose 
any problems numerically, and we typically set our initial 
conditions at z — 1000. 



6 CONCLUSIONS 

We have shown that both spherically symmetric collapse and 
cosmology can be studied in a single global coordinate sys- 
tem, which we have called the Newtonian gauge. This gauge 
considerably simplifies the task of modelling cluster forma- 
tion and its effects on the CMB. Observers can be placed in a 
region of the model corresponding to an FRW Universe, and 
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Figure 4. Density evolution for an m = 3 model. The initial density profile contains a small, compensated perturbation around the 
background level. Initially, the density falls with the expansion of the Universe. After t = 897 the streamlines turn round, the central 
region starts to collapse and the density there increases. The perturba tion evolves to give the final density distribution shown at t = 1772. 
This closely matches the observed density profiles in galaxy clusters (Dabrowski et al. 1997). The units are arbitrary. 



both the lensing and CMB temperature d ecrement are eas- 



ily mod elled. In an accompanying paper ( Dabrowski et al, 
199?|)~dletailed computations are performed with this model 
and compared with other work and with observations. Nu- 
merical integration along photon paths are performed to 
compute the lensing effects of the cluster. The effects on 
the CMB power spectrum are also computed by evolving a 
simulated background (with Gaussian fluctuations) through 
the cluster to the observer. 

A study of Swiss Cheese models within the Newtonian 
gauge revealed the surprising result that the models only 
yield unique predictions because of the particular way the 
initial data is given. A new formula for the CMB tempera- 
ture decrement in a Swiss Cheese model was presented, and 
it was argued that the number of arbitrary parameters in 
such models limits their usefulness. 

The problems with Swiss Cheese models are overcome 
in a simplified model in which a small density and veloc- 
ity perturbation is introduced at recombination and allowed 
to grow to give a compensated density profile. The centre 
of this density distribution models a cluster with a density 
profile close to that of the King profile. Photons can then 
be sent through this density distribution to model the ef- 
fects of clusters at arbitrary redshifts. The lack of pressure 



support means that the model probably over-estimates the 
perturbations induced in the CMB. This will certainly be 
the case for nearby clusters which are virialised and do not 
have large infall velocities. Recent work on structure forma- 
tion, however, has suggested that clusters with large infall 
velocities could well have formed around redshifts of 2 and 
higher. If so, the model presented in this paper could prove 
very useful for calculating the effects of such clusters. 
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From this we can extract T/r, and from all of the above 
measurements we can reconstruct F, r and u. 

In the homogeneous region, however, this final step is 
not possible because the right-hand side of ( |A3| ) vanishes. 
This prevents one from determining r from local measure- 
ments alone, which is reassuring because in the homogeneous 
region all local observables should be functions of time only. 



APPENDIX A: LOCAL OBSERVABLES 

It is instructive to see how one might measure F, r and it 
in the collapsing region with only local measurements of the 
gravitational field, and to see how the method fails outside 
the cluster. The key variable is the strength of the radial 
tidal force. For two objects a (small) proper distance d apart 
the relative acceleration of one from the other is 

a = d (inp - ^) . (Al) 

By measuring the force necessary to prevent this tidal ac- 
celeration and comparing this with the local density it is 
possible to calculate the local value of M/r 3 . Of course, in 
a homogeneous region this would be just 4np/3. 

Next we measure the variation of M/r 3 with time for 
freely-falling observers. This goes as 

D_(M\__3Mu 

Dt \r 3 J ~ r 3 r' ( ' 

so measurements of this lead to a determination of u/r. 
It follows from the definition of T (j^ that the value of 
(r 2 — l)/r 2 can also now be determined. Finally, we need 
to measure the local variation of M/r 3 with proper distance 
on a constant timeslice. This goes as 
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